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Abstract 

We present a bilinear Hirota representation of the N = 2 supersymmetric extension 
of the Korteweg-de Vries equation. This representation is deduced using binary Bell 
polynomials, hierarchies and fermionic limits. We, also, propose a new approach for the 
generalisation of the Hirota bilinear formalism in the N = 2 supersymmetric context. 


1 Introduction 

The study of exact solutions of completely integrable supersymmetric systems is of current 
interest in modern mathematical physics research. In particular, the N = 2 supersymmetric 
extension of the Korteweg-de Vries (KdV) equation [I] has been largely studied in terms of in- 
tegrability conditions, exact solutions and symmetry group structures Ejiiiiiiiisiiiiisiisiiin]- 
The equation is described by a bosonic superfield A defined on the superspace 2 [U of 
local coordinates (x,t, 0i,02)- The variables {x,t) are the usual Euclidean space-time coor¬ 
dinates and (01,02) are real Grassmann coordinates satisfying the usual anti-commutation 
relations 

0102 -b 0201 = 0, 02=6»|=O. (1) 

The superfield A satisfies the Labelle-Mathieu N = 2 supersymmetric extension of the KdV 
equation [1] 


At — (,~Axx -b (a -b 2)ADiD2A -b (a — l)(DiA)(^D2A) -b aA^)x, (2) 

where a S R is a real parameter, Di, D 2 are the covariant derivatives defined as 

D,=dg,+eA, i = l,2 (3) 

and satisfy Df = = dx- The bosonic superfield A can be decomposed m using a Taylor 

expansion as 


A(x,t,0i,02) = u{x,t) + 01^l{x,t) + d 2 ^ 2 {x,t) - dld 2 v{x,t), (4) 

where m, v are complex-valued even functions and q, ^2 are complex-valued odd functions. 

Labelle and Mathieu [T] showed that equation (I2|) is completely integrable for the special 
choices a = —2,1,4. This fact suggests that, for these special values, the supersymmetric 
KdV equation possesses travelling wave and multi-soliton solutions Him 0 0 [7]. An 
algebraic direct method to find such solutions is described by the Hirota bilinear formalism. 
This method was used numerous time for non-supersymmetric integrable evolution equation 
to construct soliton and similarity solutions, Backlund and Darboux transformations and 
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to obtain integrability conditions [T2]. Carstea has adapted this formalism to fV = 1 su¬ 
persymmetric extensions m such as the KdV, modified KdV and Sine-Gordon equations. 
The generalisation of this formalism to N = 2 extensions has been confronted to numerous 
difficulties. Zhang and al. [5] used the strategy of decomposing equation into two N = 1 
equations for which the Hirota formalism is well adapted. The way of achieving this is to 
re-write the bosonic superfield A given in (|3]) as 

A{x, t, Oi,02) = Ao{x, t,0i) + 02^{x, t, 0i), (5) 

where Aq{x, t,0i) = u{x, t) + 0iii{x^ t) and S(x, t,0i) = ^ 2 ( 2 ;, t) + 0iv{x, t) are, respectively, 
bosonic and fermionic superhelds of {x,t,0i) € In order to use the bilinear Hirota 

formalism m, we have to re-write the superhelds Aq and in terms of dimensionless 
bosonic superhelds. This is done by dimensional analysis following the fact that equation 
([2]) is invariant under the dilatation vector held [2] 

xdx + itdt + ]^0ide^ + ]j^02de^ - AOa- ( 6 ) 

This vector held shows that, under the transformation 

{xA,0u02,A) (Aa;,A3t,A^0i,A^d2,A-iH) , (7) 

equation (ED is invariant, where A is a free parameter. Under these transformations, we 
deduce the dimension of these quantities 

[H]=-l, [5,] = -l, [Di] = -\ (8) 

and this allows us to re-write the superhelds Hq and S as 

Aq = S = DiPj;, (9) 

where B = B{x,t,0i) and p{x,t,0i) are bosonic superhelds such that [B] = [p] = 0. Intro¬ 
ducing the superheld A, given by ([SD together with the relations (HD, in equation ([SD, we get 
two N = 1 supersymmetric equations 

Bt = —Bxxx-t-{ 0 , + 2)BxPxx + {<1 — ^){DiBx){DiPx) + aB^^ (10) 

Dipt = -Dipx XX + iPxxDiPx - (a + 2 )BxDiBxx + (1 - a)BxxDiBx + SaBlDipJ^ll) 


The study of these two equations is the main object of the paper for the special cases 
a = —2,1,4, i.e. for the cases for which equation ([JD is completely integrable [I]. Indeed, 
we will give a bilinear Hirota representation of these equations using different approaches 
using the binary Bell polynomials |14j . hierarchies [HlElllo] and fermionic limits [4]. The 
later approach is simple, it relies on taking = ^2 = 0 in the representation of the bosonic 
superheld A given in (j4D. The binary Bell polynomials have, recently, found a connection 
with the Hirota bilinear formalism [14] and this connection will be used throughout this 
paper. 

The one-variable Bell polynomials Y are dehned as 




( 12 ) 


where are integer constants and / = /(a;, t, 0i) is a bosonic superheld. Using the polyno¬ 
mials Y, we dehne the binary Bell polynomials y as 


ykxX,ktt,kidi (2C1, m2) 


Ykx.x,ktt,ki6i {fka 


xkftki 9 




(13) 


where the different derivatives of / are replaced by the superhelds mi and m 2 following the 
procedure 


Wi,kxxk,thei if b i® 

W 2 :kxxhthei if kx + h + ki is even. 


(14) 
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Note here that we are using the notation ~ The link with the 

Hirota bilinear formalism is given by 

ykMMe. = In .W 2 = Hfg)^ = (/ • g), (15) 

where the Hirota derivative is defined as 






■' {dt - t, 6»i)g(a;', t', 


( 16 ) 

These relations will be used to transform equations (Hoj and (HU) into binary Bell poly¬ 
nomial equations. To achieve this, we will need auxiliary tools such as hierarchies of the 
supersymmetric KdV equation ([2]) and fermionic limits. 

This paper is divided as follows. In the following three sections, we give a bilinear rep¬ 
resentation of the supersymmetric KdV equation ([2]) for, respectively, a = 1, a = 4 and 
a = — 2. In section 2, we use directly the binary Bell polynomials to get a general Hirota 
formulation. In the a = 4 case, we use the binary Bell polynomials and the Two-Boson su¬ 
persymmetric equation laii], which are members of the same hierarchy, to obtain a Hirota 
equation. In section 4, we get the bilinear representation using fermionic limits and retrieve 
the bosonic Miura transformation [12] linking a solution of the KdV equation with the mod¬ 
ified KdV equation. The last section addresses the open problem of generalizing the Hirota 
bilinear formalism to the N = 2 supersymmetric context. We exhibit this generalisation 
threw the supersymmetric KdV equation with a = 1. We conclude the paper with some 
future outlooks and remarks. 

The novelty of this paper is based on the use of the binary Bell polynomials (fT^ to get 
a bilinear representation of the supersymmetric KdV equation ([2|). As of today, the Hirota 
formulation of the supersymmetric KdV equation with a = —2 was an open problem and 
here we propose a partial answer to this question. We also, for the first time, generalize the 
Hirota bilinear formulation to iV = 2 extensions of certain integrable systems. 


2 The supersymmetric KdV equation with a = 1 

In this section, we directly use the binary Bell polynomials to get a bilinear representation 
of the supersymmetric KdV equation ([^ with a = 1. In this case, the two iV = I equations 
(HU and (HU reduces to 

Bt = —Bxxx + ‘iBxPxx + B^, ( 17 ) 

Dipt = -Dip^ XX + ipxxDip,-W,DiB,, + WlD,p,. ( 18 ) 

We can notice from these equations that the superfield B is associated to odd numbers of 

derivatives while the superfield p to even numbers. This fact is compatible with the binary 
Bell polynomials [14] . Indeed, we have the expressions 

yt{cB,dp) = cBt, y 3 x{cB,dp) = cB^^^ + 3cdB^p^:c + (19) 

from which we easily deduce an equivalent representation of equation HU given as 

y(^H,-p)+V3x(^B,-p) = 0. (20) 

For the second equation HU, we have the following binary Bell polynomials, assuming 
equation HU is satisfied, 

{iB, -p) = -Dipt + B^x^DiB - 3Bj;p^^DiB - BIDiB, (21) 

yxxx0i {iB, -p) = -HiPxxx + BIDiB + 3HxPxxT>i.B - -H ZBlDip^, 

+ 3p XX Dipx - 3BxDiB XX 
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( 22 ) 


and it is direct to show that equation (TTSl) is equivalent to the binary Bell polynomials 
equation 

yte^ {iB, -p) + ya:xx9i {iB, -p) = 0. (23) 

Using the link between the binary Bell polynomials and the Hirota derivative (fT5l) , we obtain, 
casting the change of variables iB = \n(f /g) and —p = ln(/g), the bilinear representation 
of the supersymmetric KdV equation with 0 = 10 given as 

{Vt+Vl)if-g) = 0, S,{Vt+Vl){f-g) = 0. (24) 


3 The supersymmetric KdV equation with a = 4 

In this case, we give a Bell polynomial perspective to the supersymmetric KdV equation 
with 0 = 4 using its integrable hierarchy [SI [in]- Unlike the a = 1 case, the Bell polynomial 
approach may not be directly applied to equation (nni and m with 0 = 4. Indeed, these 
equations are explicitly given as 

Bt = —Bxxx + ^BxPxx + 3{DiBx){DiPx) + 4:B^, (25) 

Dipt = -Dip^ XX + 3p 

XX ^iPx ^^X^l^XX 'iBxxDiBx + 12B^DiPx, (26) 

and the terms {DiBx){Dipj;), B^xDiBx are incompatible with the definition of the binary 
Bell polynomials. These terms are incompatible in the sense that the superfield B is asso¬ 
ciated to odd number of derivatives in the binary Bell polynomial perspective. So, we have 
to find an other way of writing these problematic terms. One way of achieving this is by 
considering the integrable hierarchy associated to the supersymmetric KdV equation with 
0 = 4. One member of this hierarchy is the Two-Boson supersymmetric equation [^ |S] for 
which its flow commutes with the flow associated to equation ([2]) with o = 4. Regarding 
this matter, the Two-Boson supersymmetric equation is given by 

Wt2 = i.-Wx + w'^+2Dip)^, (27) 

Pt2 = {px + 2.wp),:, (28) 

where w and p are, respectively, bosonic and fermionic superfields. Making use of the 
dilatation invariant vector field, we cast the change of variables w = Cx and p = Digx, 
where C and g are dimensionless bosonic superhelds. In this case, after integration with 
respect to the variable x, equations (E71) and (1^ reduce to 

Ct2 = —Cxx + + ‘2'qxx^ (29) 

Diqt2 = DiQxx + 2,CxDiqx. (30) 

At first sight, this system is again incompatible with the binary Bell polynomials, but, taking 
m = 2q — C and n = C, these equations are given as 

nt2 = mxx + nl, (31) 

Dimt2 = DiUxx+ 2.nxDimx- (32) 

We thus observe that m is associated to even numbers of derivatives, while n to odd numbers. 
From the binary Bell polynomials 

yt2in,m) = nt2, yxxeAn.m) = nixxDin + nlDin + DiUxx + 2.nxDimx (33) 
yxx{n,m) = mxx + nl, yt2eAn,m) = Dimt2 + mxxDin + nlDin, (34) 

it is direct to see that the system of equations (l3T|) and (l32|) is equivalent to the system of 
binary Bell polynomials 

Vt 2 (n,m) - Vxx(?T-,m) = 0, VtsSi («-, ™) - (n, m) = 0. (35) 
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Making the change of variables n = ln(//( 7 ), m = ln(/g) and using the link between the 
Bell polynomials and the Hirota derivative (USD, we get the bilinear representation of the 
Two-Boson equation given as 

{Vt,-Vl){f-g) = 0, S^{Vt,-Vl){f-g)=0. (36) 

Let us now focus our attention on the first equation (|25|) of the supersymmetric KdV equation 
for a = 4. We have, taking n = aB and m = /3p, the following binary Bell polynomial 
expression: 


yt2x{n,m) = yt2x{o:B,Pp) = aB^xx + iafiBxPxx + a^Bl + 2af5{DiBx)(DiPx), (37) 

where we have retrieve the previously problematic term {DiBx){DiPx). Supposing that 
equation (1251) has the binary Bell polynomial representation 

—yt{aB,Pp) = 'jyxxxiaB, Pp) + 5yt^x[aB, Pp), (38) 

a 

for a, /3, 7 , (5 S C constants, we directly find that 

z 3z 

a = 2i, P = -2, 7 = g, ^ s' 


and, thus, equation (ESI) as the following Hirota bilinear formulation, taking 2iB = ln(//g) 
and -2p = ln(/g), 

{Vt + -Vl + \VxV,)j {f-g)= 0. (40) 

The Bell polynomial analysis of the second equation (E51) is similar as for the first one. 
Indeed, it lies on the following binary Bell polynomial expressions: 


-yte^{aB,Pp) 


^yxt2eiiaB,Pp) 
(^-^i Pp) 


= Dipt - 2BxxxDiB + 12BxPxxDiB + 6{DiBx){Dipx)DiB 

+ SBlDiB, 

= ‘2‘BxxxDiB — 8{DiBx){DiPx)DiB — 12BxPxxDiB — SB^DiB 
T DiPxxx T ^BxxDiBx “t“ GBxDiBxx ^‘^B^DiPx ‘^PxxDiPx^ 
= 2BxxxDiB - 12BxPxxDiB - SB^DiB + DiPxxx - QpxxDiPx 
+ QBxDiBxx — ^‘^B^DiPx- 


From these expressions, it is elementary algebra to show that equation (1261) is equivalent to 
the binary Bell polynomials equation 


yte^{2iB, -2p) -I- \yxxx0t{‘^iB, -2p) -|- -yxt2eA‘^iB, -p) = 0 


(41) 


and, using the same change of variables for B and p as for equation (ESI) , we get the bilinear 
representation 

5i [vt + \vl + \VxVt)j {f-g) = 0. (42) 

Hence, equations (l40l) and (1^ represent the Hirota formulation of the supersymmetric KdV 
equation with a = 4 [5], where the functions B, p, f and g are related as 



p = 




(43) 
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4 The supersymmetric KdV equation with a = —2 


As of today, a Hirota bilinear representation of the supersymmetric KdV equation with 
a = —2 was an open problem [^. In this section, we propose a partial answer to this 
problem by considering the fermionic limit [3] of the superfield A given in Q. This means 
that we take = ^2 = 0 in (|4]) and we consider bosonic superfield A of the form 


A{x,t, 81 , 02 ) = u(x,t) - 6 i 02 v(x,t) (44) 

solution of equation @ with a = — 2. In this case, the bosonic complex-valued functions u 
and V satisfy the system of partial differential equations 

Ut = y^xxx bu Ux, (45) 

Vt = -Vxxx + &VVx + QUxUxx - du^Vx - l2uUxV. (46) 


We can make some observations on this system: equation (HSl) is the bosonic modified KdV 
equation for which a Hirota formulation is known and, as a second observation, we have 
that, taking u = 0, this system reduces to the bosonic KdV equation [12]. Again, at first 
sight, this system is incompatible with the definition of the binary Bell polynomials. To 
solve this problem, we use a Miura-type transformation relating the supersymmetric KdV 
equation with a = — 2 to the equation m 

Qt = —Qxxx + + -^{D2Qx)iD2Q)Qx — -^{D2Qx){DiQ){DiD2Q) 

- \{D 2 Q){D^Qx){D^D 2 Q) + \{DiQx){DiQ)Qx, (47) 

where Q is a bosonic superfield. The equation (l47l) is the first non-trivial flow of a iV = 2 
supersymmetric hierarchy, as shown by Tian and Liu [10] . and the Miura-type transformation 
is explicitly given as 

A =—D 1 D 2 Q-\-^{D 2 Q){DiQ). (48) 

Using the Taylor expansion Q{x,t, 81 , 02 ) = qo{x,t) — 8 i 62 qi 2 {x,t) (fermionic limit), the 
Miura-type transformation (1481) is equivalent, in components, to 


1 


u = -qi2, 


— ^(912 + Oo , x ) 


1 


(49) 


Note that setting u = 0 in these transformations leads to the classical Miura transformation 
V = \qQ^x ~ \oo,xx relating a solution of the modified KdV equation to a solution of the KdV 
equation m- This observation is compatible with the fact that the bosonic complex-valued 
functions go and gi 2 satisfy the decoupled system of partial differential equations 


qo,t 

912,t 


qo,xxx T 


1 


3 

Oo,x, 


qi2,xxx 2 ^ 12912 ,: 


(50) 

(51) 


These two equations are all of the modified KdV-type and we know that they possess a 
Hirota bilinear representation [12] . Indeed, they are equivalent to the binary Bell polynomial 
equations 


yt(l/2qo,/3qo) +yxxx(l/2qo,/3qo) - 3yx(l/2qo,/3qo)yxx(i/2qo,/3qo) = 0, (52) 
yt(i/2pi2,7Pi2) + yxxx('i/2pi2,7Pi2) - 3yx(i/2pi2,7Pi2)yxx(i/2pi2,7Pi2) = 0, (53) 


where gi 2 = Pi 2 ,x, P and 7 are arbitrary constants and go and pi 2 are auxiliary and arbitrary 
bosonic functions. Using the change of variables 



/3go = ln(/g). 



7Pi2 = ln(/5). 


(54) 
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we get the Hirota bilinear representation 


{Vt+Vl){f-g) = 3AP,(/-5), Vl{f-g) = Xfg, (55) 

iVt+VDCf-g) = iXD^Cf-g). Vl{f ■~g) = \f~g, (56) 

where A and A are arbitrary parameters and the functions u and v take the explicit forms 

u =In , V = {d^\nf/gf - {d^\nf/gf - dllnf/g. (57) 

In the reminder of this section, we give plots of different soliton solutions. For a 1-soliton 
profile, the functions /i, gi, fi and gi, solution of the bilinear equations (I55|) and (l56l) with 
A = A = 0, may be chosen as 

/i = l + e’’, gi = l-e’', /i = l + ie’', gi = l-ie'^, (58) 

where rj = kx — K^t and f] = kx — k^t. In the case of a 2-soliton profile, the functions / 2 , 52 , 
/2 and 52 given as 

/2 = 1 + + 6 "= + ^ 126 ''!+’'=, 32 = 1 - , (59) 

/2 = 1 + - ii 2 e^i+^", 32 = 1 - (60) 

solves the bilinear equations (l55l) and (1551) for A = A = 0, where gi = KiX — n^t, fji = kiX — k^t 
for i = 1, 2, Ai 2 = ( j and A 12 = ( . In the figures below, we use the notation 





V(m,n) = {dx ^-0-fmlgmf “ {dx ^^fnlgnf “ 9^ ln/m/ 3 m (61) 


and we have made the choices /t = l,K=|,/ti = |, K 2 = ^, Ki = f and K 2 = §■ Figure I 
represents the functions U(i,i) and W( 2 , 2 ); figure 2 gives the behavior of the functions U(i_i) 
and V( 2 . 2 ), while figure 3 plots the functions U(i_ 2 ) and U( 2 q). 



Figure I: The plot of the functions U(i,i) and M( 2 , 2 )- The dashed curve correspond to time 
t = — 2 , the thick curve to time f = 0 and the other curve to time t = 2. 


5 The N = 2 formalism: a first approach 

In this section, we make an attempt to formalize the Hirota bilinear approach in the N = 2 
supersymmetric context using the Bell polynomials. We will illustrate this new approach 
on the supersymmetric KdV equation ([2]) with a = I. Using the change of variable A = Mx 
and integrating once with respect to the variable x, the supersymmetric KdV equation for 
a = I reads as 

Mt = -Mxxx + iMxDiD2Mx + Ml, (62) 
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Figure 2: The plot of the functions U(i i) and V( 2 , 2 )- The dashed curve correspond to time 
t = —2, the thick curve to time t = 0 and the other curve to time t = 2. 




Figure 3: The plot of the functions ^(i 2 ) and V( 2 ,i)- The dashed curve correspond to time 
t = —2, the thick curve to time t = 0 and the other curve to time t = 2. 


where M = M{x, t, 0i, ^ 2 ) is a dimensionless bosonic superfield. In order to adapt the Hirota 
bilinear formalism in the N = 2 context, we make the following observation 

=d^^D, = d-^de,+e, = dl{D-^D-^M,) = C,,, (63) 

where C = and d~^ is the integration operator with respect to x. In this 

setting, the equation (I62[) can be re-written as 

Mt = -M,,, + 3T4a, + m3 (64) 

and we observe that the Bell polynomial approach can be directly used. Indeed, the above 
equation is equivalent to the binary Bell polynomial equation 


yt{iM, -C) + J,,,(iM, -C) = 0 (65) 

which, using the identification iM = \n(f / g) and —C = hi{fg), can be written as 

{Vt+Vl)if-g) = 0 , ( 66 ) 

where / = f{x,t, 61,62) and g = g{x, 1,61,62) are bosonic superfields. The identification 
imposes a further constraint given as 

D2Di\n{fg) = (67) 

and it can be shown that this relation is equivalent to the bilinear Hirota equations 


S 2 Si{f ■ f) = 2iff,^, S 2 Si{g ■ g) =-2igg,^. (68) 


In order to study this additional constraint, we consider the following Taylor expansions 

ln(/) =/o + 0 i/i + 6*2/2 + 0i02/i2, ln(g) = go + 6 *i5i + 6*252 + 6*16*2512 ( 69 ) 















and, once introduce in (I67L we get 


lii(/) — /o + (^1 ~ *^2)/i + *^1^2/0,XI In(ff) — 90 + (^1 + *^2)<?i ~ 'i(^i& 29 o,x- ( 70 ) 

For the 1-soliton solution, we can choose 

^ _j_ Q^g77+0iCl+^2C2+^1^27ni2 g = \ -\- ^g'^ + 01^'1+02^'2+^1^2'71'12 (71) 

as a solution of the bilinear equation (l66l) for a, 6 ,mi 2 ,ni 2 bosonic quantities, Cii C 21 
V 2 fermionic quantities and rj = kx — K^t. These free parameters have to be determined 
in order that the constraint (p71) be satisfied. This can be done by substituting the above 
expressions for / and g in (1701) . We get the explicit forms 

fo = ln(l + ae’'), /i = aCi—- C 2 =-iCi, mi 2 = in, (72) 

1 + ae^ 

go = ln(l + 6e’'), gi = bvi , V 2 = ivi, ni 2 = -in (73) 

1 + be^ 

and, as a consequence, the components of the superfield A defined in (|3]) are given as 

u = i{go,x-fo,x), = i{gi,x - fi,x), 6 =-(/i.x+ gi.x), v =-(/o,xx + go.xx)- (74) 

In this section, we have produced, for the first time, a N = 2 supersymmetric Hirota 
representation of the KdV equation with a = 1. Indeed, the Hirota formulation is described 
by equations (1551) and (155)) . We have thus partially solved the open problem of finding a 
supersymmetric N = 2 generalisation of the Hirota bilinear formalism [5] . 

6 Future outlooks and remarks 

In this paper, we have presented a systematic way of finding the Hirota bilinear represen¬ 
tation of the N = 2 supersymmetric KdV equation using its decomposition into two N = 1 
supersymmetric equations and the Bell polynomials. For the completely integrable cases 
a = 1,4, we have obtain a complete representation using the integrable hierarchy and the 
binary Bell polynomials. For the a = —2 case, it was an open problem to find an Hirota 
formulation. Here, we have proposed a partial answer to this question using fermionic limits 
and have retrieved the well known Miura transformation relating the KdV and modified 
KdV equations. It still remains to find its general representation. 

An other open problem was to find a N = 2 generalisation of the Hirota bilinear formal¬ 
ism. We have, for the first time, succeeded in given such a representation for the supersym¬ 
metric KdV equation with a = I. The main idea in this construction was to re-write the 
operator D 1 D 2 as a second derivative with respect to x of a given quantity. This had the 
effect of transforming the N = 2 equation into a ’’new” equation involving only derivatives of 
the bosonic variables (x, t) for which the binary Bell polynomials could be directly applied. 
This as led to its iV = 2 Hirota representation. 

Our future goals is to generalize the results of section 5. This new approach avoids 
treating a N = 2 equation as two iV = I supersymmetric equations and this allows us to 
obtain a N = 2 bilinear representation of the equation. As a final example to illustrate the 
efficiency of this new proposed procedure, we consider the N = 2 supersymmetric extension 
of the potential Burgers equation 

Mt = D 1 D 2 MX + ^Ml (75) 

As in section 5, we re-write the quantity D 1 D 2 MX as Cxx where C = D^^D^^Mx and, as 
a consequence, the Burgers equations [3] reads as 

Mt = Cxx + ImI (76) 
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This equation may be cast into a binary Bell polynomials equation as 


2Xyt{XM, 2X^C) - y^4XM, 2X^C) = 0 , (77) 

where A is a free complex parameter. Making the change of variables AM = ln{f/g) and 
2X^C = \n{fg), we get the Hirota bilinear equation 

{2XV,-Vl){f-g) = 0 (78) 


together with the constraint 

D 2 D 1 ln(/g) = 2A5, In , (79) 

which can be re-write into a Hirota bilinear form as 

5251 (/■/) = 4 A// 3 ;, S 2 Si{g ■ g) =-‘iXgg:^. (80) 

In conclusion, equations (l76l) and (1801) constitute the Hirota bilinear representation of the 
supersymmetric potential Burgers equation (1761) . 
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